Nonextensive Pesin identity. Exact renormalization group analytical results for the 
dynamics at the edge of chaos of the logistic map. 
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We show that the dynamical and entropic properties at the chaos threshold of the logistic map 
are naturally linked through the nonextensive expressions for the sensitivity to initial conditions and 
for the entropy. We corroborate analytically, with the use of the Feigenbaum renormalization group 
(RG) transformation, the equality between the generalized Lyapunov coefficient \ q and the rate of 
entropy production K q given by the nonextensive statistical mechanics. Our results advocate the 
validity of the g-generalized Pesin identity at critical points of one-dimensional nonlinear dissipative 
maps. 
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The nonextensive generalization of the Boltzmann 
Gibbs (BG) statistical mechanics [1], [2] has recently 
raised much interest and provoked considerable debate 
[3] as to whether there is firm evidence for its applicabil- 
ity in circumstances where a system is out of the range of 
validity of the canonical BG theory. Recognition and un- 
derstanding of the existence of such a limit of validity is a 
major concern in the development of present-day statis- 
tical physics. So, increased attention has been drawn to 
the examination of physical situations that do not satisfy 
the customary BG equilibrium conditions, e.g. insuffi- 
cient randomness and limited or nonuniform motion over 
pertinent phase space that result in anomalous dynamical 
properties [4] , [5] . Within these several types of physical 
phenomena and connected model system properties [2], 
the critical states of one-dimensional non-linear dissipa- 
tive maps stand out. It has become apparent [2] that 
they are bona fide examples where the predictions of the 
nonextensive generalization of the BG statistical mechan- 
ics are appropriate. Here we give analytical proof and nu- 
merical corroboration of the until now only conjectured 
[6] (/-generalization of the Pesin identity between entropy 
production rate and Lyapunov exponent at the onset of 
chaos of unimodal maps. Beyond its intrinsic value as a 
means to study the dynamics of incipient chaotic states, 
this result has important implications regarding the suit- 
ability of the Tsallis entropy in describing a specific phys- 
ical situation unreachable by BG statistics. 

Recently, the predictions of the nonextensive theory 
have been rigorously proved for the pitchfork and tan- 
gent bifurcations and for the edge of chaos of logistic-type 
maps by means of the analytic renormalization group 
(RG) derivation [7] of the q-exponential expression, 
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for the sensitivity to initial conditions £t containing the 
entropic index q and the ^-generalized Lyapunov coeffi- 
cient X q . Eq. (1) has been proposed [6] as the nonex- 
tensive counterpart to the usual exponential sensitivity 
£t = exp(Ait) to initial conditions which prevails when 
the ordinary Lyapunov coefficient Ai is nonvanishing. 
(The BG exponential form is recovered when q — > 1). 

Pioneering work [8], [9] on the dynamics at the edge 
of chaos of the logistic map was directed at the deter- 
mination of the fluctuation spectrum of the algebraic 
Lyapunov coefficients X q . Here we focus on the entropic 
properties of trajectories at this state with the idea of in- 
vestigating the existence of a generalized Pesin identity. 

As a significant windfall of our RG calculations we now 
know the expressions for X q at the mentioned critical 
states of logistic-like maps [7], These expressions have 
been interpreted in terms of the fixed-point map param- 
eters and corroborated numerically via a priori calcula- 
tions [7] . Specifically, for the edge of chaos /ioo of the lo- 
gistic map Xq (and q) are simply given by X q = In a/ In 2 
(and q = 1 — In 2/ In a) where a is the Feigenbaum's 
universal constant that measures the power-law period- 
doubling spreading of iterate positions. Having reached 
this level of knowledge on £ t and X q it is only natural 
to enquire about its relationship with the entropic prop- 
erties of trajectories at /x^. The Pesin formula that re- 
lates the Kolmogorov-Sinai (KS) entropy K-i (described 
below) and the Lyapunov coefficients of nonlinear maps 
has become an extremely useful tool for the quantitative 
analysis of the dynamics of chaotic states. This formula 
embodies the all-important connection between the loss 
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of information measured by ICi and the Lyapunov coef- 
ficients A^ for chaotic states [10], [11]. The general in- 
equality ICi <J2 ^ where the sum is over the a£° > 
reduces for one-dimensional systems to the Pesin identity 
id = Ai, Ai > 0. 

So, as a starting point we consider the ^-generalized 
rate of entropy production K q , defined via K q t = S q (t) — 
S q (0), t large, where 



S q = ^Pt^q (-7 



(2) 



is the Tsallis entropy, and where ln 9 y = (y 1 q — 1)/(1 — 
q) is the inverse of exp q (y). In the limit q — > 1 
K q becomes K x = t^S^t) - 5i(0)] where Sx(t) = 

- Sill Pi (*) ln In E q- (2) Pi(t) is the probability 
distribution obtained from the relative frequencies with 
which the positions of an ensemble of trajectories occur 
within cells i = 1,...,W at iteration time t. The initial 
conditions for these trajectories have a prescribed distri- 
bution Pi(0) and the phase space into which the map is 
defined is partitioned into a large number W of disjoint 
cells of sizes k. As a difference from K\ the KS entropy 
has a more elaborate definition since it considers the en- 
tire trajectories from their initial positions to the time 
limit t — > oo [11]. The relationship between K,\ and K\ 
has been investigated for several chaotic maps [12] and 
it has been established that the equality K\ = K\ occurs 
during an intermediate stage in the evolution of the en- 
tropy Si(t), after an initial transient dependent on the 
initial distribution and before an asymptotic approach to 
a constant equilibrium value. The g-generalized KS en- 
tropy JC q is defined in the same manner as /Ci but with 
the use of Eq. (2). Here we look into the analogous in- 
termediate regime in which K, q — K q . As we shall see 
below it turns out to be sufficient to evaluate the rate 
K q for uniform initial distributions defined in a parti- 
tion of equal-sized cells to establish the validity of the 
conjectured [6] form K q = X q of the Pesin identity. 

Next we recall that the logistic map, f^(x) = 1 — 
fix 2 , —1<x<1, exhibits several types of infinite se- 
quences of critical points (with Ai = 0) as the control 
parameter fi varies across the interval < \i < 2. One 
such important sequence corresponds to the pitchfork bi- 
furcations [10]. The accumulation point of the pitch- 
fork bifurcations is the Feigenbaum attractor that marks 
the dividing state between periodic and chaotic orbits, at 
/loo = 1.40115.... A measure of the amplitudes of the pe- 
riodic orbits is defined by the diameters d n (n = 0, 1, ...) 
of the 'bifurcation forks' at the 'super-stable' periodic or- 
bits of lengths 2" that contain the point x = 0. These 
super-stable orbits occur at ~p n < and approach ^ 
as ~fi n — /ioo ~ 5~ n (n large) where 6 = 0.46692... is one of 
the two Feigenbaum's universal constants. The diameter 
d n = fjj- (0) is the iterate position closest to x = in 



such 2"-cycle. For large n these distances have constant 
ratios d n /d n+ i = —a, where a = 2.50290... is the 2nd of 
the Feigenbaum's constants [10]. For clarity we use only 
the absolute values of positions so below d n means \d n \. 

The main points in the following analysis are: 1) We 
determine the evolution of all orbits at /Uqo, i.e. those 
with initial conditions Xi„ belonging to the attractor, to 
the repeller, and to all other positions. 2) We obtain £ t 
for any x in and find the remarkable property that Eq. 
(1) holds in general with the same fixed values for q and 
X q up to a time T — 2 N where N ~ — \nx in / In a. 3) We 
observe that the position-independent form found for X q 
implies that ensembles of trajectories expand in such a 
way that a uniform distribution of initial conditions re- 
mains uniform for all later times t < T where T marks the 
crossover to an asymptotic regime. 4) As a consequence 
of this we establish the identity of the rate of entropy 
production K q with X q . We corroborate numerically all 
our findings. 

In Fig. 1 we show the absolute values of the positions 
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of two trajectories of the logistic map at 

^oc in logarithmic scales. One corresponds to x in = 0; 
and the other one to Xi n ~ 0.56023..., close to a repeller, 
the unstable solution of x = 1 — ji^x 2 . The 1st tra- 
jectory maps out the attractor while the 2nd exhibits a 
long-lived transient stretch as it falls into it. In Ref. [7] 
it was shown that a distinct fraction of the positions of 
the trajectory with Xi n = consists of subsequences gen- 
erated by the time subsequences Tk — 2™ + 2"~ fc , with 
k = 0, 1, ... and n > k, and each of these exhibits the same 
power-law decay. The main subsequence (k = 0) can be 
expressed (via a time variable shift to = To — 1) a s the 
q-exponential x to — expg(Agto) with Q = 1 + ln2/lna 
and Aq = — In a/ In 2. The positions in the subsequences 
can be obtained from those belonging to the supercycles 
at ~p n < yUoo. In particular, the subsequence for k = 
was identified as the diameter sequence = d n = 
a~ n . This property was shown to imply Eq. (1) with 
q = 1 — In 2/ In a and X q = In a/ In 2. Notice that 
q = 2 - Q as cxp q (y) = 1/ exp Q (-y). 

We can extend considerably the above results. First, 
we note that the entire attractor can be decomposed into 
position subsequences associated with the same power- 
law decay, and that all subsequences are generated by 
the time subsequences = (2k + l)2"~ fc . (The first po- 
sition at r = 2k + 1 of each of the first 8 subsequences 
can be identified among those labelled in Fig. 1). We 
make use of this time-position classification to point out 
that the positions x Tk of all trajectories with xi n < a~ n \ 
m = n — k, are given by 
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where we have neglected terms of 0(a 3m xf n ) and where 
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gk = g( 2k+1 ^(x) is the (2k + l)-th composition of the 
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fixed-point map g(x). Both g(x) and gk{x) are solutions 
of the RG doubling transformation consisting of func- 
tional composition and rcscaling, R/(x) = af(f(x/a). 
Eq. (3) is obtained by keeping the first two terms in 
the power-series expansion of gk(x) [10] followed by use 
of gk(x) — a m g^ \x/a m ) and the change of variable 
Considering a pair of initial conditions j/j„ 
and Xi n in Eq. (3) yields 

x Tk (yin) ~ x Tk (x m ) = [x 2k +i(y m ) - x 2 k+i{x m )]a m (4) 

For each subsequence k, the sensitivity £ tfc , defined as 

\ x t k (y in) Xt k {Xin)\ 



I»»-Xi»h0 \x tk =o(yin) -X tk =o(x in )\ 



(5) 



can be written, with use of the shifted time variable 
t k = T~k — 2k — 1 (n > fc), and observing that a m — 
(1 +t fe /(2fc + l)) lnQ / ln2 , as the g-exponential 



tt k = cx P(? \fh k 



(6) 



where q = 1 - In 2/ In a and A^ fc) = lna/((2fc + l)ln2). 
The crossover time T = 2 N is determined from the con- 
dition Xi n < a ~ N m Eq. (3). 

Let us consider next an ensemble of M trajectories with 
initial positions xi n uniformly distributed along the in- 
terval [1 — 1,1], for transparency with 1 — I < g^(0). 
An arbitrary partition of [1 — 1, 1] is made with a cer- 
tain number I of nonintersecting intervals of lengths li, 
i = 1,2,..., I, with I = J2i^i- For I sufficiently small, 
under Tk = (2k + l)2 n ~ k iterations the lengths /, trans- 
form, according to Eq. (4), as — a m li. Since we also 
have 1 {t ^ = a m l, we observe that the interval ratios re- 
main constant, that is k/l = Vf k ' /l^ Tk \ Thus, the initial 
number of trajectories within each interval Mlijl remains 
fixed for all times r < T, with the consequence that the 
original distribution is uniform for all times t <T. 

We can now calculate the rate of entropy production. 
This is more easily done with the use of a partition of 
W equal-sized cells of length I. Fig. 2 provides a strik- 
ing corroboration of the time constancy of uniformity. 
This figure shows in logarithmic scales the evolution of a 
distribution pi (r) of positions of an ensemble of trajecto- 
ries at Hoc beginning from a uniform distribution Pi(l) of 
initial positions contained within a single cell of size I ad- 
jacent to x = 1. If we denote by W tk the number of cells 
that the ensemble occupies at the shifted time tk and 
by Ax tfe the total length of the interval these adjacent 
cells form, we have W tk = Ax tk /l and in the limit I — > 
(since W tk = (Ax tk / Ax tk =o)(Ax tk=0 /l)) we obtain the 
remarkably simple result W tk — £t k ■ As the distribution 
is uniform, and recalling Eq. (6) for £ tfc , the entropy is 
given by S q (t k ) = In, W tk = Xq k h k , while 



as Wt k =o = 1. Eq. (7) is our main result. The numerical 
results shown in Fig. 3 substantiate and bring to light in 
a dramatic manner the validity of the g-generalized Pesin 
identity at /Zoo. 

An interesting observation about the structure of the 
nonextensive formalism is that the cquiprobability en- 
tropy expression ln g Wt can be obtained not only from 
S q in Eq. (2) but also from 
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where Sq = S 2 -q = S q . The inverse property of 
the g-exponential reads \n q y = — ln 2 - 9 (l/y) for the q- 
logarithm and as pointed out introduces a pair of con- 
jugate indexes Q = 2 — q with the consequence that 
while some theoretical features are equally expressed by 
both S q and Sq some others appear only via the use 

of either S q or Sq. For instance, the canonical ensem- 
ble maximization of Sq with the customary constraints 

Y^iLi Pi — 1 an d Y^iLi Pi e i = U, where ej and U are con- 
figurational and average energies, respectively, leads to a 
Q-exponential weight (with Q > 1 when q < 1). On the 
other hand the partition function is obtained via the op- 
timization of S q . The mutual Eqs. (2) and (8) elegantly 
generalize the BG entropy. 

We summarize our arguments and findings. Critical 
states with vanishing Ai in dissipative one-dimensional 
nonlinear maps display power law £t. It is natural to 
expect this to imply a corresponding power law rate of 
entropy production linked to the dynamics of ensembles 
of trajectories. A connection between these two proper- 
ties suggests an extension of the Pesin identity JCi = Ai, 
Ai > that incorporates the case Ai = 0. But, interest- 
ingly, to study this situation formally one is required to 
develop a theory beyond the usual BG scheme for chaotic 
states that supplies generalizations for both the KS en- 
tropy and the Lyapunov exponent. One known source is 
the nonextensive statistics constructed around the Tsallis 
entropy S q as this offers specific and practicable expres- 
sions for these quantities. To make a meaningful analysis 
of this problem it is indispensable to carry out an explicit 
a priori determination of all quantities involved and here 
we have proved that this is indeed the case for a specific 
but prototypical example, the onset of chaos of the lo- 
gistic map. We have shown [7] that the Feigenbaum RG 
method, from which the static fixed-point solution g(x) 
was originally obtained, is also capable of delivering dy- 
namical properties, most visibly the sensitivity £ t . It is 
important to stress that the derivation of £ t does not use 
in any way the nonextensive formalism and for this reason 
it constitutes an independent corroboration of the expres- 
sion for £ t suggested by this theory. With an analytical 
expression for X q in hand a parallel expression for the rate 



3 



of energy production K q was here obtained from two in- 
gredients: i) a distribution function Pi(t) of positions for 
an ensemble of trajectories and ii) the Tsallis entropy S q . 
Our main result, the generalized identity K q — X q , neces- 
sitates that the equiprobability entropy has the precise 
analytical form \n q W t (with q = 1 — In 2/ In a) and to 
this extent distinguishes the Tsallis expression S q from 
other alternatives, including the BG Si. 

We have shown for the first time that the Pesin identity 
holds rigorously, albeit in a generalized form, for incip- 
ient chaotic states. Because the entropic index q (as is 
the case of X q and its identity K q ) is obtainable in terms 
of the Feigenbaum's a we are able to address the much- 
asked question regarding the manner in which the index 
q and related quantities are determined in a physical ap- 
plication. The generic chaotic state is that associated 
to Ai > 0, but it is evident that the critical state with 
Ai = carries with it completely different physics. The 
analysis was specifically carried out for the Fcigcnbaum 
attractor of the logistic map but our findings clearly have 
a universal validity for the entire class of unimodal maps 
and its generalization to other degrees of nonlinearity. 
In a more general context our results indicate a limit of 
validity to the BG theory based on Si and the appro- 
priateness of the nonextensive S q for this kind of critical 
dynamic states. 
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FIG. 1. Absolute value of two trajectories at /loo in log- 
arithmic scales. Empty circles correspond to Xi n = (the 
numbers label time t = 1, 16). Small dots correspond to 
x in — 0.56023..., close to a repeller, the unstable solution of 
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FIG. 2. Time evolution, in logarithmic scales, of a distribu- 
tion pj (r) of trajectories at fj,^ . Initial positions are contained 
within a cell adjacent to x — 1 and j counts the consecutive 
location of the occupied cells at time r. Iteration time is 
shown for the first two subsequences (k = 0, 1). 



FIG. 3. Numerical corroboration (full circles) of the gen- 
eralized Pesin identity — at fioo. On the verti- 
cal axis we plot the q- logarithm of £ tfc (equal to \q h h) and 
in the horizontal axis S q (equal to Kgt). In both cases 
q = 1 — In 2/ In a = 0.2445.... The dashed line is a linear 
fit. In the inset the full lines are from the analytical result 
Eq. (6). 
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